function od=optical_density(trans_signal,in_signal)
od=-log(trans_signal./in_signal);


end